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Abstract 

The tunneling conductance spectra of a normal metal / insulator / singlet 
superconductor is calculated from the reflection amplitudes using the Blonder- 
Tinkham-Klapwijk (BTK) formulation. The pairing symmetry of the super- 
conductor is assumed to be d x 2_ y 2 + is, or d x 2_ y 2 + id xy . It is found that in 
the (d x 2_ y 2 + is)-wave case there is a well defined conductance peak in the 
conductance spectra, in the amplitude of the secondary s-wave component. 
In the (d x 2_ y 2 + id xy )-w&ve case the tunneling conductance has residual val- 
ues within the gap, due to the formation of bound states. The bound state 
energies depend on the angle of the incident quasiparticles, and also on the 
boundary orientation. On the basis of this observation an electron focusing 
experiment is proposed to probe the [d x 2_ y 2 + id xy )-w&ve state. 
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I. INTRODUCTION 



Two decades ago, Blonder, et. al. used the Bogoliubov-de Gennes (BdG) equations 
to calculate the tunneling conductance of normal metal / s-wave superconductor contacts, 
with a barrier of arbitrary strength between them, in terms of the probability amplitudes 
of Andreev || and normal reflection. In the Andreev reflection process an electron inci- 
dent, in the barrier can be reflected as an electron (normal reflection), reflected as a hole 
without changing its momentum (Andreev reflection), it can also be transmitted into the 
superconductor as an electron-like, hole-like quasiparticle. 

Recently the BTK theory was extended by several groups to consider the anisotropy 
of the pair potential. In <i-wave superconductors the pair potential changes sign under a 
90°-rotation. So under appropriate orientation of the a-axis of ci-wave superconductor the 
transmitted quasiparticle feel different sign of the pair potential. This results in the forma- 
tion of bound states within the energy gap, which are detected as peaks in the conductance 
spectra. In <i-wave superconductor the peak exists at E = for a great range of angles 
of incidence of the incoming electron. This range depends on the surface orientation. |3j 
In particular for (110) surfaces the peak exists at E = for all angles of incidence, and 
disappears for the (010) or (100) surface. 

In the presence of another barrier inside the normal metal additional subgap bound states 
exist due to multiple Andreev reflections @||. The same phenomenon occurs in <i-wave 
superconductor / insulator / d-wave superconductor. In these systems the quasiparticle 
current has been examined by several groups using BTK formalism with recursive 

relations for the determination of the probability amplitudes. 

There is a competition between different pairing symmetries in the bulk. The coexistence 
of a sub dominant order parameter in bulk depends on the strength of the secondary order 
parameter attractive interaction relative to the attractive interaction in the dominant pairing 
channel. When the secondary order parameter is strong enough a second phase transition 
occurs at a temperature T cl < T c which depends on the strength of the secondary order 



2 



parameter. Numerical results show that when such coexistence is realized the relative phase 
of the order parameters is 7r/2 leading to d x 2_ y 2 + is or d x 2_ y 2 + id xy pairing state in the 
bulk. The temperature dependence of the various thermodynamic quantities and transport 
properties change from power lows to exponential below T c \. |9|,[lIJ When the secondary 
order parameter is not strong enough, only the d^^-wave order parameter appears in the 
bulk. For (110) surfaces the d x 2_ y 2-wave order parameter changes sign under reflection at 
the surface and vanishes at the surface. On the other hand the s or d xy does not change 
sign and are not effected by the presence of the surface, so there is the possibility of their 
presence near the surface even when their attractive interaction is not strong enough for 
them to exist in the bulk |TT|JT^| . 

The presence of the secondary order parameter near a surface is manifested in tunneling 
spectra as a splitting of the zero energy conductance peak (ZEP) at low temperatures at 
zero external field and further non-linear splitting with increasing external field. |13j The 
field dependence of the splitting of the ZEP in the tunneling spectra of YBCO has been 
examined flT3| , |r5|] . The observation is consistent with a d x 2_ y 2 + is surface order parameter 
or a d x 2_ y 2 + id xy order parameter. 

In this paper we extend the BTK formula to calculate the tunneling conductance in a 
normal metal / insulator / [d x 2_ y 2 + is)-wave, or [d x 2_ y 2 + id xy )-wave superconductor. In 
particular we find that in the d x 2_ y 2 + is-state the conductance peak remains rigid at the 



energy of the subdominant (s) order parameter. [|16|,|T7| Besides in the d x 2_ y 2 + id xy state, 
there is a plateau region inside the gap due to the formation of bound states at discreet values 
of the quasiparticle trajectory angle 8, for all junction orientations. Also the evolution of 
the tunneling conductance with temperature depends on the nature of the subdominant 
order parameter. These features can be used to distinguish between states with broken 
time-reversal symmetry. 
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II. THE MODEL FOR THE NS INTERFACE 



We consider the normal metal / insulator / superconductor junction shown in Fig. |l|. 
We choose the y direction to be parallel to the interface, and the x direction to be normal 
to the interface. The insulator is modeled by a delta function, located at x = 0, of the form 
VS(x). The temperature is fixed to OK. 

The motion of quasiparticles in inhomogeneous superconductors is described by the BdG 
equations 



(1) 



H e (r)u(r)+ Jdr'A(s,x)v(r') = Eu{r) 
fdr'A*(s,x)u(r')-H* e (r)v(r) = Ev(r) ' 

where the single-particle Hamiltonian is given by Ti. e (r) = —h 2 \j 2 r /2m e + V(r) — Ep, E 
is the energy measured from the Fermi energy Ep. A(s,x) is the pair potential, after a 
transformation from the position coordinates r, r' to the center of mass coordinate x = 
(r + r')/2 and the relative vector s = r — r' . After Fourier transformation the pair potential 
depends on the related wave vector k and x. In the weak coupling limit k is fixed on the 
Fermi surface (\k\ = hp), and only its direction 9 is a variable. Also we neglect any spatial 
variation near the interface, e.g. the pair potential does not depend on x. The pair potential 
has the form: 



A fa;, 



(2) 



0, x < 
A(0), x > ' 

where 9 is the angle of the quasiparticle trajectory measured from x-axis. When a beam of 
electrons is incident from the normal metal to the insulator, with an angle 9, the general 
solution of Eqs. (Jl]), is the two component wave function, which for x < is written 
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while for x > 0, the solution is 
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where a, 6, are the amplitudes for Andreev and normal reflection, and c, d are the amplitudes 
for transmission into the superconductor as electron-like and hole-like quasiparticles respec- 
tively. In the following we assume that q e « w fc e « fc^ ~ fop. The latter approximation 
is valid within the BCS weak coupling theory. The BCS coherence factors are given by 



u. 



i = ll + jE*-\A ± (6)\yE]/2, 



(5) 



and 



[1-V£ 2 -|A ± (0)| 2 /£]A 



(6) 



The internal phase coming from the energy gap is given by <p± = [A±(9)/\A±(9)\], where 
A_)_(6 I ) = A(#) (A_(0) = A(7r — 9)), is the pair potential experienced by the transmitted 
electron-like (hole-like) quasiparticle respectively. Using the matching conditions of the wave 
function at x = 0, tf/(0) = tfjj(0) and W n {0) - %(0) = (2mV/?i 2 )^ I (0), the magnitude of 



the Andreev and normal reflection R a = \a\ and Rj, = \b\ , are obtained as [ITJ 
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(1 — Ujv)|1 — n + n_(fi_(j) 



* 12 



(8) 



|1 + (o~n - l)n + n_0_0+| 2 ' 

where n± = v±/u±. The tunneling conductance, normalized by that in the normal state is 
given by Q 



a{E) 



r_l) 2 d9a s {E^ 



J-i%d9a 



(9) 



N 



according to the BTK formula the conductance of the junction, a s (E,9), is expressed in 
terms of the probability amplitudes a, and b: W S (E, 9) = 1 + R a — R b . The transparency of 
the junction ctn is connected to the barrier height V by the relation 



where Z = 2mV/f?kp, denotes the strength of the barrier. In the Z = 0(large ctjv) limit the 
interface is regarded as a weak link, showing metallic behavior while for large Z[a^ = 0) 
values the interface becomes insulating. 
We consider the following cases: 

a) In case of d x a_ y 2-wave superconductor 

A(0) = A 1 (r)cos[2(0 -(3)1 (11) 

where (3 denotes the angle between the normal to the interface and the x-axis of the crys- 
tal. The temperature dependence of the gap follows the usual BCS relation Ai(T) = 
A ds J\ — T/T d , where T d is the transition temperature. 

b) In the [d x i_ y i + is)-wave case 

A(0) = Ai(T) cos[2(0 - p)] + zA 2 (T), (12) 

where A 2 (T) = A s \Jl — T/T s , and T s is the transition temperature for the s-wave compo- 
nent. 

c) In the (d x 2_ y 2 + id xy )-wave case 

A(0) = Ai(T) cos[2(0 - /?)] + iA 2 (T) sin[2(0 - /?)], (13) 

where the angular form of the secondary component is obtained by the substitution of (3 in 
the d x 2_ y 2-wave order parameter by (3 + 7r/4. A 2 (T) = A dxy ^Jl — T/T dxy , follows the BCS 
relation, and T dxy is the transition temperature for the d xy -wave component. 

III. TUNNELING CONDUCTANCE CHARACTERISTICS 

In Figs. 2-4 we plot the tunneling conductance cr(E) as a function of E/Aq for various 
values of Z, for different orientations (a) (3 = 0, (b) 7r/8, (c) 7r/4. The pairing symmetry 
of the superconductor is d x 2_ y 2-w&ve, with A^ = 0.7A , in Fig. ||], [d x 2_ y 2 + is)-wave, with 
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A d = 0.7A , A s = 0.3A , in Fig. |), (d x 2_ y2 + id xy )-wave, with A d = 0.7A , A dxy = 0.3A 
in Fig. |j. It is clear from these figures that the peaks are narrowed by the increase of Z. In 
this section the temperature is fixed to OK. 

For (3 = e.g. when the lobes of the dominant o?-wave component point towards the 
junction interface the position of the conductance peak, is near the energy gap A^, in all 
the above pairing symmetries. This peak is mainly effected from the bulk density of states. 

For (3 7^ another peak exists in the conductance spectra, for the d x 2_ y 2-wave, [d x 2_ y 2 + 
is)-wave cases, but its physical origin is different than that found near A^. For the <i-wave 
case this peak exists at E — for all the non zero values of /3, due to the different sign 
of the pair potential that the transmitted quasiparticles feel. However, the height of the 
conductance peak (ZEH) depends on the orientation angle j3. For a given angle j3 the ZEH 
is proportional to the range of 9 angles for which sign change occurs. This is seen in Fig. 
(c) for (3 = it/ 4 where the ZEH is maximum since for this orientation the transmitted 
quasiparticles feel different sign of the pair potential for all angles —ir/2 < 9 < n/2. On the 
other hand for f3 = ir/8 in Fig. |2| (b) the range of angles is reduced and the ZEH takes a 
lower value. 

For the [d x 2_ y 2 +is)-wave case in Fig. ^ the position of the conductance peak is shifted to 
the energy E = A s , for all values of /3. For each value of (3, its height depends on the range 
of 9 angles where the transmitted quasiparticles feel different sign of the pair potential. For 
(3 = 7r/4 the conductance peak as seen in Fig. |3| (c) at E = A s , has its maximum value since 
the transmitted quasiparticles feel the sign change of the pairing potential for all angles 9. 
This range is reduced for other orientations and for (3 = it goes to zero, as we can see 
in Fig. [3] (a). Also a subgap opens within the conductance spectra due to the imaginary 
s-wave component. Within the subgap in the tunneling limit, the tunneling conductance 
is zero, <j(E) = while in the metallic limit (Z = 0), <y(E) = 2 independently from the 
orientation, as in the s-wave case. In the Z = case the normal reflection coefficient is 
zero while the Andreev reflection coefficient is unit. In this case the charge transport into 
the superconductor is twice as large as in the normal state, for energies within the subgap 
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region. 

In the d x 2_ y 2 + id xy case, seen in Fig. |] the tunneling conductance has residual values 
within the gap for all orientations (3. In particular for (3 = 0, as seen in Fig. |] (a) in the 
tunneling limit, the conductance a(E — 0), at E — has a non zero value contrary to the 
(d x 2_ y 2 + is)-wave case where it is zero. In the d x 2_ y 2 + is case for (3 = 0, there is no angle 9 
for which the transmitted quasiparticles to experience the sign change of the pair potential, 
and the tunneling conductance goes to zero. This is different for the d x 2_ y 2 + id xy where 
for f3 = the transmitted quasiparticles feel the sign difference due to the secondary order 
parameter d xy . 

Also the zero energy conductance height evolves very differently with the orientation 
of the superconductor, for the three pairing symmetries. This is seen in Fig. |5| where the 
dependence of zero-energy conductance height on (3 is plotted, for Z = 2.5, for the three 
pairing symmetries. It is seen that for the [d x 2_ y 2 +id xy )-w&ve case (dashed line), for (3 close 
to 7r/4, the height representing the plateau like feature seen in Fig. |] is enhanced. Besides 
for angles close to zero, the height of the ZEP for the (d x 2_ y 2 + id xy )-w&ve case is reduced, 
but is not zero. Note that the height for (3 = 0, remains finite even in the large Z-limit, 
while the height in the [d x 2_ y 2 + is)-wave case goes to zero. 



IV. BOUND STATE ENERGIES 

These features are explained if we calculate the energy of the midgap state, which is 
given for large Z by the value in which the denominator of Eq. ([7|,[8]) vanishes. The equation 



giving the energy peak level is written as [16 



<f)^<f)ln + n_\ E=Ep = 1.0. (14) 

In the d x 2_ y 2-wave case, for a given angle (3, this equation has solution E — 0, for a fi- 
nite range of angles 9. For (3 = 7r/4 the solution is E = for —ir/2 < 9 < n/2, since 
n + n_\E=o — — 1, and also the transmitted quasiparticles feel different sign of the pair poten- 
tial i.e. 0-01|b=o = — 1- In the [d x 2_ y 2 + zs)-wave case for /3 — 7r/4, the solution is E = A s 
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in the 9 interval [0 : vr/2]. In this case the n + ,n_ and the internal phases are varied in a 
way that Eq. ^ is satisfied for E = A s and a midgap state is formed. When a midgap state 
exists the tunneling conductance a s (E, 9) is equal to 2 for all 9 and the peek in cr(E) seen in 
Fig. is due to the normal state conductance cn in Eq. [H], which depends inversely on 
the Z 2 for large Z. For intermediate angles ft, the peak height of the tunneling conductance 
ct(E) is proportional to 8/3, for < ft < tt/4, || and for ft = the range of 9 angles for 
which Eq. [14| has solutions collapses to zero in both symmetry states, and no bound states 
are formed. Then u{E) goes to zero as 1/Z 2 and there is no conductance peak. For energies 
different than the bound state energy E p , for large Z, W S (E, 9) is inversely proportional to Z 2 
as the ctat is and the tunneling conductance has a constant value as we can see in Fig. for 
E > 0. In the [d x 2_ y 2 + id xy )-w&ve case for fixed ft the solutions of [L4| depend both on E, and 
9, as seen in Fig. ||, where the bound state energy E p is plotted for ft = 0, 7r/16, tt/8, 7r/4, 
as a function of 9. In this case the midgap state for a given ft is formed for a pair of angles 
9, for energies within the gap. This observation can be used to explain the residual values 
of the tunneling conductance within the gap, seen in Fig. |] as follows. When a bound 
state is formed the conductance a s (E,9) is equal to 2 exactly at the bound state energy 
for the two discreet values of 9 and the peak in the cr(E) should be proportional to the Z 2 
for large Z for these values of 9. For the rest of the quasiparticle trajectory angles 9 the 
tunneling conductance <j(E), has a constant value. Thus the height for a given energy, and 
angle ft is determined from the interplay of two competitive factors, i.e., the bound state 
energy formed at a couple of 9 angles which gives a contribution proportional to Z 2 , and 
the rest of 9 angles which give a constant value contribution independent of Z. Also the 
steps in 9, in evaluating the integral in Eq. |9] are very much crucial since the calculation 
of the tunneling conductance has to be performed exactly at the bound state energy. If 
this is not the case then the peak due to the bound states in the tunneling conductance 
would have a smaller value which would depends on Z in general. We conclude that in the 
d x 2_ y 2 + id xy the discreet values of the quasiparticle trajectory angle 9, over which a bound 
state is formed, compared to the interval of 9 angles in the other two pairing states explains 
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the reduced height of the tunneling conductance within the gap. However if we calculate 
the conductance a s (E, 9), for the d x 2__ y 2 +id xy at a given j3 for a value of 9 for which bound 
state occurs then the conductance should develop a peak at the bound state energy, where 
a s (E,9) is equal to 2. For the rest of the energies, a s (E,9), goes to zero as 1/Z 2 . This is 
seen in Fig. [7] where the conductance a s (E,9) for Z = 2.5 is plotted for fixed (3 = it/A as 
a function of the energy E, for different values of the angle 9 = n/4, 3n/8, n/2 for which 
bound state occurs. We see that for 9 = 7r/4 the peak is at E — 0. However as we change 
the angle 9 towards tt/2 the peak level moves from E = to E = . 

The occurence of residual density of states in the (d x 2_ y 2 + id xy )-waye case is unaffected 



by the calculation of <j(E) including the self-consistency of the order parameter. [18| In this 
calculation an enhancement appears at E — A^ , for (3 = it/ 4. In our calculation we also 
observed a similar enhancement at E = Aj when the definition 



d9a s (E,9) cos 9 

f-ir/2 d9(JN COS 9 



°m = J "^/77r'^7" . ( 15 ) 



was used for the calculation of the tunneling conductance. The cos 9 factor was included in 
the integration formula to calculate the x-component of the tunneling spectra. Within this 
definition the bound state at 9 = contributes more (due to the cos 9 factor) that the bound 
state at 9 close to 7r/4. As seen in figure ^| the bound state at 9 = corresponds to energy 
E = Ad xy causing the peak in u{E) at E = A^. Also in a self consistent calculation the 
bound state at {6 = 7r/4, E = 0) contribute less in u{E) that that at {6 = 7r/4, E — A^) due 
to the depletion near the interface. In any case the peak near E = in (d x 2_ y 2 + id xy )- 
wave pairing state is much more suppressed that that at E = A s in (d x 2_ y 2 + id)- wave 
state. 

The angular dependence of the bound state energy for fixed boundary orientation at 
the xy plane can be used to identify the [d x 2_ y 2 + id xy )-waMe pairing state. The method we 
propose here is the two point spectroscopy described by Benistant et. al. |0|. They measured 



the reflected hole distribution along the boundary y-direction, when electrons are injected 
with certain distribution P (</>), through a point contact, at y — 0, into a normal metal of 
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thickness d attached to an s-wave superconductor. The presence of a magnetic field parallel 
to the z-axis deflects the trajectories of the electrons and leads to an asymmetric distribution 
of angles of incidence in the normal metal / superconductor interface. Also the magnetic 
field focuses the reflected holes into a second point contact which acts as a hole collector. 
Moving the second point contact around the first one or using several point contacts along 
the direction parallel to the interface we are able to measure the intensity of the Andreev 
reflected holes as a function of the y direction. In the s-wave case one observes a single peak 
called 'focusing peak' aiy — yo from the injection point at y = 0, since the Andreev reflected 
probability amplitudes are independent of the injection angle. For the (d x 2_ y 2 + id xy )-wave 
case the bound state energies, for which the reflection coefficient is equal to one, occurs 
at angles 8\ < 0, 82 > 0, for a given boundary orientation and large barrier strength Z. 
These bound states will give rise to a second peak in the hole distribution, at a different 
position, besides the one due to the focusing. The presence of the magnetic field leads to 
an asymmetric distribution of angles of incidence in the interface and the trajectory which 
corresponds to bound state at 8\ has shorter path from that at 82 and the corresponding 
injected electrons have smaller angle 0. If the angular distribution probability P(4>) of the 
injected particles is peaked at small injection angles, this will lead to larger contribution 
to the secondary peak from the bound state at 8\ then that at 82 This new peak would be 
observed for all energies for which bound states exist in a (d x 2_ y 2+id xy )-w&ve superconductor. 
In the case of a o^-^-wave superconductor, and a (d x 2_ y 2 + is)-wave superconductor the 
resonance exist only for E = 0, E = A s correspondingly. Experiments of this kind require 
high quality normal conducting crystal and point contacts for the electron injection. Any 
voltage drop has to occur at the point contact for the electrons to move ballistically in 
the normal metal. Similar procedure has been proposed from Honerkamp and Sigrist [2(| to 
discriminate between unitary and nonunitary triplet states for the superconductor Sr 2 Ru04. 
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V. TEMPERATURE DEPENDENCE OF THE TUNNELING SPECTRA 



At finite temperatures the tunneling conductance is calculated from the relation [^T 

/oo 
dE 
-oo 



df(E + eV) 



dE 



j:i 2 /2 d6a s (E, 



/-tt/2 



(16) 



eV is the electron energy and f(E) is the Fermi function f(E) = l/(er E + 1), (3 = l/fc^T. 
In the case of a two component order parameter, we assume that below a surface transition 
temperature a subdominant order parameter can develop which breaks spontaneously the 
time reversal symmetry. Its amplitude is bellow the value for the formation of spontaneously 
broken time reversal symmetry state in the bulk. [|l2j The temperature dependence of the 



pair potential amplitude is assumed to obey the usual BCS relation. As a consequence 
under the coexistence of the secondary component the T c 's for the dominant c?-wave Td and 
the subdominant s(d xy ) components, T s (Td xy ) directly correspond to the amplitude of the 
attractive interaction in each case. 

Fig. |8] shows the tunneling conductance a(eV) for different temperatures T /Td = 
0.1, 0.2, 0.3, 0.4, in the large barrier strength limit Z = 10, (3 — tt/4. The pairing symmetry 
of the superconductor is d^.^-wave in Fig. |8| (a), (d x 2_ y 2 + is)-wave, with T s = 0.37^, in 
Fig. H (b), (d x 2_ y 2 + id xy )-w&ve, with Td xy = 0.3Td in Fig. § (c). It is seen that due to the 
thermal occupation of states contributing to the tunneling current, the peaks are getting 
broadened as the temperature increases. In the 6^2 _ y 2 -wave case, as seen in Fig. |8] (a) 
the ZEP is suppressed when the temperature increases and disappears almost at the criti- 
cal temperature. This feature of the calculated spectra is consistent with the experimental 
results of YBa2Cu307_5 observed by low-temperature scanning tunneling spectroscopy. (§] 
The evolution of the conductance spectra with temperature is qualitatively similar with the 



calculation including the self-consistency. |22| On the other hand in the (d x 2_ y 2 + is)-wave, 



seen in Fig. |8] (b) the tiny subgap of the order of A s = 0.3Ao at T = disappears with the 
increase of the temperature. For T > T s it follows the usual d x 2_ y 2 -wave like dependence. 
In the (d x 2_ y 2 + id xy )-w&ve, shown in Fig. § (c) the zero energy height is suppressed with 
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the increase of the temperature. For T > Td xy the temperature dependence of the spectra 
for the (d x 2_ y 2 + id xy )-wa.ve state is similar to the d x a_ y a-wscve. 

In Fig |9] we plot the ZEH as a function of temperature for the three pairing states. For 
the d x 2_ y 2-wave case the ZEH behaves as T _1 . For the (d x 2_ y 2 + zs)-wave case the ZEH 
increases up to T = 0.2T^ and then decreases with increasing the temperature. For T > T s 
it follows the d x 2_ y 2-wave behavior. The downturn of the ZEH at low temperatures in the 
(d x 2_ y 2 + is)-wave case corroborates with the ZEP splitting and has also been observed 
experimentally (see Fig.l in Ref. |13|). In the (d x 2_ y 2 + id xy )-wave case the ZEH decreases 
with T in different scales for T < Td xy and T > indicating the different pairing states. In 
all cases the transition from the d x 2_ y 2-wave to the [d x 2_ y 2 + is)-wwe or [d x 2_ y 2 + id xy )-w&ve 
is continuous. 

In the metallic limit [Z = 0) (not presented in the figure) the tunneling conductance at 
eV = decreases, as the temperature increases, from its zero temperature value a(eV) = 2, 
to the normal state value er(eV) = 1 at the transition temperature. The variation with 
T for the d x 2_ y 2 + is (d x 2_ y 2 + id xy )-w&ve for T < T s (Td ) deviates from the d^.^-wave 
behavior. In both cases where time reversal symmetry is broken, a change of slop exists in 
the a(eV = 0) vs T diagram, at the subdominant order parameter transition temperature. 
However in this case the variation with T is similar for d x 2_ y 2 + is, d x 2_ y 2 + id xy and thus it 
can not be used to discriminate between the two pairing states. 

VI. CONCLUSIONS 

We calculated the tunneling conductance in normal metal / insulator / anisotropic su- 
perconductor, using the BTK formalism. We showed that the conductance peak for (110) 
surface orientation, in a d x 2_ y 2-wave superconductor, appears in zero energy, and is shifted 
according to the amplitude of the secondary order parameter in the [d x 2_ y 2 +is)-wave case. 
In the (d x 2_ y 2 + id xy )-wave case the tunneling conductance has residual states within the 
energy gap. These are due to the formation of bound states at discreet values of the tra- 
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jectory angle 9 for each boundary orientation angle (3, for energies within the gap. These 
bound states explain both the residual states within the subgap and also the small height 
of the conductance within the subgap region. The calculation of the conductance a s (E,9) 
for given boundary orientation at an incident angle 9 for which bound state occurs shows 
an enhancement at the bound state energy. The energy dependence of the bound state on 
9 can be used within the method of electron focusing to detect the [d x 2_ y 2 + id xy )-wscve 
state. In such a case besides the focusing peak another peak exists in the reflected hole 
distribution spectrum for all energies of the injected electrons less than the amplitude of 
the secondary order parameter. This peak should also be observed for the cL^^-wave and 
(d x 2_ y 2 + is)-wave cases but only at the energy E = 0, E = A s respectively. 

The zero energy conductance peak decreases as T _1 with increasing the temperature and 
disappears almost at the transition temperature for the G^.^-wave case. The temperature 
dependence of the ZEH deviates from the usual T _1 behavior of the d x 2_ y 2, in the case 
where a subdominant surface order parameter is developed, for T < T cl , where T cl is the 
transition temperature for the subdominant order parameter. These features can be used 
to distinguish between time-reversal broken symmetry states. 

Throughout this paper the spatial variation of the dominant order parameter near the 
surface which depends on the boundary orientation is ignored for simplicity. As a conse- 
quence, since the nucleation of the secondary order parameter near the surface depends on 
the strength of the dominant one, the spatial variation of the secondary order parameter is 
also ignored. We expect more drastic changes when the orientation is (3 = tt/A where the 
suppression of the dominant order parameter is more significant. However, since the features 
presented here are intrinsic and are generated by the existence of surface bound states, the 
essential results do not change qualitatively. 

Also we assumed perfectly flat interfaces in the clean limit, so any impurity scattering 
and the effect of the surface roughness are ignored. Generally surface roughness will lead to a 
statistical distribution of the outgoing trajectories, and will alter the results presented. The 
effect of surface roughness on the tunneling effect in interfaces between normal metals and 
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superconductors with time reversal symmetry broken, has been studied. |23| It is found that 
in the d x 2_ y 2 + id xy case additional bound states are formed due to the surface roughness. 
Also in d x 2_ y 2-wave superconductor, the ZEP may appear even for (100) interfaces with 
surface roughness [fT^ ]. 

Also in a more realistic treatment of the problem, one has to take into account also the 
thickness of the barrier. In that case additional resonances are expected in the tunneling 
spectra due to multiple Andreev reflections within the barrier, besides the ones due to the 
bound states. 
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FIG. 1. The geometry of the normal metal / insulator / superconductor interface. The vertical 
line along the y-axis represents the insulator. The arrows illustrate the transmition and reflection 
processes at the interface. In this figure, (3 is the angle between the normal to the interface and 
the a-axis of superconductor, and 6 is the angle of the incident electron beam and the normal. On 
the top, the ci-wave order parameter is shown. 
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FIG. 2. Normalized tunneling conductance <j{E) as a function of E/Aq for Z = (solid line), 
Z = 2.5 (dotted line), Z = 10 (dashed line), for different orientations (a) (3 = 0, (b) ir/8, (c) 7r/4. 



The pairing symmetry of the superconductor is d x 2_ y 2, = 0.7Aq. The temperature is T = 0. 
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FIG. 3. The same as in Fig. 2. The pairing symmetry of the superconductor is d x 2_ y 2 + is, 
A d = 0.7A , A s = 0.3A . 
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FIG. 4. The same as in Fig. 2. The pairing symmetry of the superconductor is d x 2_ y 2 + id xy , 



A d = 0.7A , A dxy = 0.3A . 
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FIG. 5. Normalized tunneling conductance a for E = as a function of (5 for Z = 2.5, 
and T = 0. The pairing symmetry of the superconductor is d x i_ y i (solid line), A^ = 0.7Ao, 
d x 2_ y 2 + is (dotted line), = 0.7Ao, A s = 0.3Ao, and d x 2_ y 2 + id xy (dashed line), A^ = 0.7Ao, 
A dxy = 0.3A . 
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FIG. 6. Bound state energy E p , for T = 0, versus the quasiparticle angle 9 for different 
orientations (3 = 0, vr/16, tt/8, tt/4. The pairing symmetry of the superconductor is d x 2_ y 2 + id xy 
with A d = 0.7A , A dxy = 0.3A . 
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FIG. 7. Conductance a s (E,6) for Z = 2.5, and T = 0, as a function of i? for fixed angle 
/? = 7r/4, at different angles 6 = vr/4, 3ir/8, tt/2, for which a bound state is formed at a different 
value of E. The pairing symmetry of the superconductor is d x 2_ y 2 + id xy with = 0.7Aq, 



A 



0.3A . 
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FIG. 8. Normalized tunneling conductance a versus the applied voltage eV, for different 
temperatures T/T d = 0.1, 0.2, 0.3, 0.4. The barrier strength is Z = 10, and the junction orientation 
is fixed to (5 = tt/4. The pairing symmetry of the superconductor is d x 2_ y 2 in a), d x 2_ y 2 + is with 
T s = 0.3T d in b), and d x 2_ y 2 + id xy with T dxy = 0.3T d in c). 
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FIG. 9. Normalized tunneling conductance a for eV = as a function of the temperature 
T /Td for Z = 10, and (3 = 7r/4. The pairing symmetry of the superconductor is d x 2_ y 2 (solid line), 
d x 2_ y 2 + is with T s = 0.3Trf (dotted line), and d x 2_ y 2 + id xy with Td xy = 0.3Trf (dashed line). 



25 



